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Recent measurements suggest the possibility to exploit ionic liquids (ILs) as smart lubricants for 
nano-contacts, tuning their tribological and rheological properties by charging the sliding interfaces. 

Following our earlier theoretical study of charging effects on nanoscale confinement and squeezout 
of a model IL, we present here molecular dynamics simulations of the frictional and lubrication 
properties of that model under charging conditions. First we describe the case when two equally 
charged plates slide while being held together to a confinement distance of a few molecular layers. 

The shear sliding stress is found to rise strongly and discontinuously as the number of IL layers 
decreases stepwise. However the shear stress shows, within each given number of layers, only a weak 
dependence upon the precise value of the normal load, a result in agreement with data extracted from 
recent experiments. We subsequently describe the case of opposite charging of the sliding plates, 
and follow the shear stress when the charging is slowly and adiabatically reversed in the course of 
time, under fixed load. Despite the fixed load, the number and structure of the confined IL layers 
changes with changing charge, and that in turn drives strong friction variations. The latter involve 
first of all charging-induced freezing of the IL film, followed by a discharging-induced melting, both 
made possible by the nanoscale confinement. Another mechanism for charging-induced frictional 
changes is a shift of the plane of maximum shear from mid-film to the plate-film interface, and 
viceversa. While these occurrences and results invariably depend upon the parameters of the model 
IL and upon its specific interaction with the plates, the present study helps identifying a variety 
of possible behavior, obtained under very simple assumptions, while connecting it to an underlying 
equilibrium thermodynamics picture. 

PACS numbers: 68.35.Af,68.08.De,62.10.-l-s,62.20.Qp 


INTRODUCTION 

Ionic liquids (ILs) - organic salts that are liquid at 
room temperature - are of considerable and increasing 
physical and technological interest. ILs are characterized 
by negligible vapor pressures, high temperature stability 
and a wide electrochemical window [T]. Moreover their 
physical properties can be widely tuned by changing the 
molecular structure of the cation-anion pairs [2]. Many 
ILs strongly adhere to solid surfaces [3] and can form 
wear-protective films capable of resisting much higher 
loads than those of molecular lubricants [1]. For this 
reason among others their use as boundary lubricants 
has been recently pursued in experimental studies un¬ 
der nanoscale confinement, employing both the surface 
force apparatus (SFA) [HIS], and atomic force microscopy 
(AFM) [SfHOj. as well in in a number of computer simu¬ 
lations at various levels of idealization Our spe¬ 

cific interest focuses here on the exploration of the fric¬ 
tional changes induced by electrical charging of the plates 
which confine the IL, a subject only partially covered by 
previous theoretical works [iMn]- Besides lubrication, 
the dynamical behavior of ILs under charging is of ad¬ 
ditional interest in the field of supercapacitors [Mini- 
Experimentally, there is a number of nano-frictional stud¬ 
ies of surfaces charged with respect to a reference elec¬ 
trode [SIIIIIH]- Not so between two oppositely charged 


electrodes. Recent macroscopic friction studies have be¬ 
gun to appear for IL sliding under charging conditions 
[Bun]. Here however we shall restrict ourselves to fric¬ 
tion under conditions of nanoscale confinement. 

At the molecular level, ILs confined between hard plates 
become structured in the form of layers that run parallel 
to the plates, not unlike other liquids but further charac¬ 
terized by the charge order typical of molten salts, with 
an alternation of positive and negative ion layers, and 
an interlayer separation that corresponds to the ion pair 
size [HIIlliniEMTl. Plate charging is expected to give 
rise to rearrangements of the IL layering and ordering, at 
least near the confining plate-IL interfaces, with a conse¬ 
quent and probable change of the lubrication properties 
under plate sliding. Charging-induced friction changes 
have already been to some extent detected and described 
experimentally li M- Not unexpectedly, they are re¬ 
ported to depend on the particular kind of IL and of 
confining surfaces chosen. For example, a strong reduc¬ 
tion of the friction of an AFM silica colloid probe was 
found by negatively biasing a Au(lll) surface [8] im¬ 
mersed in [Pyi^iFAP], However, a similarly strong fric¬ 
tion drop occurred at positive bias upon sliding a sharp 
AFM tip on highly oriented pyrolytic graphite immersed 
in a [HMIm]FAP IL 0. 

This diversity of behavior, reflecting the variety of ILs 
and the different nature of the confining and sliding sur- 
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faces, suggests some flexibility in the corresponding the¬ 
oretical modeling, aimed at a broad exploration of how 
some of this variety of behavior could be addressed with 
a handful of model parameters, rather than describing in 
details one or another particular case. The minimal IL 
model and a natural first choice is a simple molten salt, 
such as liquid NaCl. Fedorov and Kornyshev [T3] did 
that, their model consisting of two spherical Lennard- 
Jones (LJ) particles with unequal radii further endowed 
with opposite charges. The frictional behavior of that 
model was recently investigated by Fajardo et al. [T3] 
who noted interesting similarities with data by Li et 
al. [9]. However, a parallel study which we conducted of 
this type of model under simple confinement and squee- 
zout m revealed that the confined liquid film evolved 
far too sharply from a layered liquid for neutral plates at 
large spacings, to a strongly crystallized rocksalt struc¬ 
ture solid under narrow confinement and/or plate charg¬ 
ing. On top of that, the simple molten salt model wet¬ 
ted the plates much too readily, even given a reasonable 
choice of parameters. Both of these aspects, excessive 
crystallization and excessive wetting tendencies, repre¬ 
sent a severe oversimplification by comparison with real 
IL lubricants whose wetting capability is generally more 
modest, and where, owing to far more complex molecu¬ 
lar structures, the layering and solidification which occurs 
under strong confinement does not usually imply cystal- 
lization, and is replaced by glass-like disordered struc¬ 
tures. That complex behavior is of course much better 
captured by fully realistic IL model simulations, of which 
there are several good examples in literature, mm but 
whose difficulty and intrinsic complexity make the studies 
much more work-intensive while at the same time empha¬ 
sizing aspects that are to some extent undesirably specific 
to each case. 

A previous study of our own IL model m was 
designed to bridge the gap between these two extremes, 
either oversimplification or overspecificity. We showed 
that a minimal modification of the charged LJ model, 
consisting of a neutral ’’tail” rigidly attached to one of 
the ions (the cation in our case, but the anion could 
have been equivalently chosen) could go a long way 
in our desired direction, making the structural and 
squeezout properties of the confined IL model consid¬ 
erably more realistic while still not specific. While the 
anion and the charged part of the cation retain a strong 
Coulomb correlation, the neutral tail is generally and 
merely pushed out of the way, resulting in poor tail-tail 
correlations. The tail-induced disorder, along with the 
steric hindrance they introduce, impede overall rocksalt 
crystallization with drastic changes in the ILs solidi¬ 
fication, also introducing new unexpected structures 
as well as the desired glassy aspects and wettability 
features. These structural features were found to evolve 
and change under plate charging, foreshadowing a 
variability that could be explored through a variation of 


the model parameters. At the same time, this IL ’’tailed 
model” (TM), to be further detailed below, is simple 
enough to permit a novel, quantitative calculation and 
characterization of the confined IL enthalpy versus 
interplate distance, and - of crucial importance here - 
its variations under plate charging. The minima of the 
enthalpy-distance curves identify stable or metastable 
layering states of the confined IL, showing a perfect fit 
with the evolution in the structural layering oscillations 
of density and charge observed in simulations of the 
confined IL under variable interplate separation and 
plate charge. The layering transitions during squeezout, 
well established in real IL experiments OEaES!, are 
reproduced by the TM model and directly related to 
jumps between enthalpy minima, while these minima 
evolve and cross values under closing of the gap between 
plates, or under variable plate charging [TS]- With that 
very useful characterization, the TM model provides a 
ready tool for a fresh study of some generic effects of 
electrical plate charging on the nanoscale sliding friction 
of confined ILs. 

In this work we first introduce and detail the TM 
model and the chosen confining geometry. We then 
analyze the sliding simulations and discuss how the 
structural changes, induced by a negative plate charging, 
modify the frictional properties. Interesting observations 
emerge examining the shear velocity profile inside the IL 
film as a function of the plate charging. First, we observe 
a switching of the shear band plane from the plate-film 
interface, typical of low friction, to the mid-film region, 
typical of larger friction. A second important effect 
is the charging-induced solidification/melting process 
described in our previous, static study, and typically 
ignored in literature. Charged plates generally increase 
the IL wetting, giving rise to ’’electrowetting”. But, as 
most ILs are not far above freezing at room temperature, 
increased wetting will in turn encourage solidification, 
the capability to withstand load, with a large effect 
on friction. Conversely, plate neutrality and partial 
dewetting facilitates melting of the IL film, leading to 
squeezout under even moderate load. 

The final and most specific part of this work is the in¬ 
vestigation of the effects on friction of opposite charging 
(equivalent to an applied interplate voltage) of the two 
confining plates. A given charge state of the plates is 
accompanied by a strong structuring of the IL near and 
between the plates. A reversal of plate charges is followed 
by an accompanying destructuring and even melting of 
that layering structure, prior to the reversal of local IL 
ordering and eventual re-solidification. The sliding fric¬ 
tion between plates undergoes a corresponding charge- 
dependent evolution as a consequence. This kind of phe¬ 
nomena, only poorly explored so far, constitute a fresh 
prediction of this study. 
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FIG. 1: (Color online) Sketch of the simulation geometry with 
open boundaries along x and a directions and periodic bound¬ 
aries along y, i.e. the direction along which the sliding force 
is applied through a spring moving with speed V, a situation 
conceptually similar to an SFA setup. The width L of the 
plates along x is 20 nm. By applying an external z-oriented 
force to the top plate, both at V =0 and at V >0, the size D 
of the gap between the plates can change since the liquid can 
flow in and out the lateral reservoirs. The shapes of the anion 
(blue) and cation (red) are sketched in the lower left corner. 


MODEL, SIMULATION DETAILS, SYSTEM 
PREPARATION 

To pursue our desired broad scenario of possible 
charging effects on the lubrication properties of nano- 
confined ILs we adopt the simple TM model m 
which deliberately neglects the fine structure of the 
ions retaining only minimal features common to most 
of the ILs, such as the anion-cation asymmetry, their 
possible amphiphilic character and their generally 
irregular shape. The TM model is a two-component 
charged Lennard-Jones (LJ) system where anions and 
cations not only have different radii but where a tail 
is attached to the cation. The cation is therefore a 
dimer consisting of a positively charged, small-sized LJ 
head, rigidly bound to an equal size, neutral particle, 
see Fig. The effect of the tail upon the ILs wetting 
and confinement properties is important, as mentioned 
in the introduction and recently shown in ref. m, 
where the details about the adopted potentials and 
the corresponding parameters are also specified and 
discussed. Essentially the tails give the molecules a 
larger gyration ratio in the liquid state, and prevent low 
temperature ordering replacing crystallization with a 
glass-like solidihcation - all features that are common 
and important in real ILs. All the molecular dynamics 
(MD) simulations were performed using the LAMMPS 


code |23]. The long range Coulomb interactions have 
been treated in reciprocal space using a particle-particle 
particle-mesh solver. Canonical ensemble configurations 
were generated by means of a Langevin thermostat 
directly applied to the lubricant molecules. The plates 
were treated as rigid bodies, the lower one fixed and 
the upper one subjected to a directed load force F/v 
and driven to slide along (x,y). Their structure consists 
of rigid close-packed triangular lattices of LJ particles 
with nearest neighbor spacing 0.52 nm, mimicking the 
periodicity of a mica surface. A graphite, or graphene 
plate, which unlike mica may be electrically charged 
has a periodicity of about 0.246 nm, which is not too 
different from half that value. In most of our simulations 
the liquid is confined between two identical plates with 
a modest surface charge density q = subse¬ 

quently increased to q = —12^C/cm'^. As a reference, 
the surface charge density quoted for mica in contact 
with H 2 O is -33.4 ^C/cw?, a much larger value [30] . 
The state of charge of a mica surface in contact with an 
IL will certainly differ, but it is still generally believed 
that some charging remains. Since nature as well as the 
algorithm used to treat the long range interaction in 
reciprocal space require the overall system to be neutral, 
when endowing both plates with a total net nonzero 
charge we correspondingly add extra ions of the opposite 
charge to compensate. 

The lateral drive is actuated through a spring k 
connected to the top plate and pulled at constant 
velocity V. The same value V = ^flrajs was adopted 
throughout unless another value is explicitly stated. 
The instantaneous frictional force opposing the motion 
is measured through the elongation of the driving 
spring as Fi(t) = k{Vt - YcM{t)), where Ycmit) is 
the y coordinate of the top plate center of mass. The 
relevant quantity in friction is eventually the shear stress 
a =< Fl > /A where A is the area of the contact. In our 
case A = 177 nm^ is a constant, so presenting < Fl > or 
a is just the same. In experiments however the contact 
area is not generally constant, so that the shear stress 
must be extracted before making comparisons. 

The open geometry described in Fig. has been chosen 
because it permits particles to be squeezed out or sucked 
in from two lateral IL droplets. These droplets thus 
serve as liquid reservoirs allowing the number of ions 
effectively confined inside the gap to adjust and change 
dynamically depending on the loading, charging, and 
sliding conditions, realizing a kind of effectively grand- 
canonical configuration. The distinguishing feature of 
the open boundary geometry employed here, as opposed 
to several previous studies which used closed periodic 
boundary conditions, is to address a situation much 
closer to real SFA or AFM experiments, permitting in 
particular squeezeout transitions and the consequent 
transverse (z) and in-plane {x, y) reordering. 



4 


Before the sliding simulations begin, we prepared equi¬ 
librium starting configurations. In our protocol, the 
starting configurations were obtained by filling the open 
gap between initially distant plates, in any desired state 
of charge, with IL, and then reducing the gap width D 
with D = 0.11 m/s, generally slow enough to enact a 
reasonably adiabatic squeezout, as was found in Ref. m- 
The adiabaticity is facilitated by our chosen temperature 
T = 225K, substantially higher the TM bulk melting 
temperature T^eit — ISORT m, but as will be discussed 
below it is not automatically guaranteed, because both 
confinement and charging favor an increase of viscosity 
and a tendency to solidify. Above a minimal interplate 
distance Dc — Anm and for small plate charge, the sim¬ 
ulated IL is completely liquid, even if structured in close 
vicinity of the plates owing to their rigid wall nature and 
also to their charging. Interestingly, this near-plate local 
structure develops not only along z, i.e., normal to the 
interface, but also along {x,y), parallel to the interface. 
Parallel IL ordering, rather more difficult to detect exper¬ 
imentally, is currently being pursued and detected with 
AFM tools [211 [ 21 ] . As H is further and gradually de¬ 
creased, the confinement effects increase, the liquid-like 
diffusion diminishes and below a critical film thickness 
Dc there is an effective interlocking of the structure ema¬ 
nating from the two facing plates, leading to a freezing of 
the IL into a solid-like layered arrangement with an odd 
number Niayer of alternating charge layers - when the 
two plates are neutral or identically charged. While no 
claim can obviously be made to represent a realistic situ¬ 
ation, we note that IL confinement between mica plates, 
typically used in SFA friction experiments, is known to 
behave precisely like that [U | 22 ] . 

As was shown in the previous paper m our choice of 
scheme and geometry permits additional thermodynam¬ 
ical insights. Treating the interplate separation D as the 
external control variable, a free enthalpy H {D, Fn, q) can 
be calculated first estimating the free energy W(D) by 
integration of the average force < F{D) > measured be¬ 
tween the plates as their distance D is reduced. 

pOO 

W{D) = < F{z) > dz (I) 

Jd 

where < F{z) > is the average z-oriented force exerted 
by the IL on the plates while they are a distance z apart. 
After that the free enthalpy at given load and charge can 
be obtained as 

HiD,FN,q) = W{D)+FND, ( 2 ) 

Free enthalpy curves are displayed in Fig[^ for two val¬ 
ues of normal load F/v, and predict the load-induced 
transition Niayer —>■ Niayer — 2 layers. The plate-plate 
free enthalpy generally displays as a function of D sev¬ 
eral odd-A^ minima whose relative values depend on plate 
charging and load. Once specified charge and load, one of 


them is the absolute minimum, therefore the thermody¬ 
namically stable layered state, and all others are higher, 
metastable states, which may still be very long lived. If 
the layered state is obtained by applying a force rather 
than by forcing a distance D, the precise number of layers 
attained in the course of a given simulation (or experi¬ 
ment) will be history-dependent, in that more than one 
observed thickness may correspond to the same applied 
force. This picture is quite close to what is observed in 
SFA and AFM [H [5] [23] ■ In addition to the alternating 
charge layer structure, we had previously found that the 
confinement-induced solid structure of the TM model, al¬ 
though far from completely cystalline, carried a second 
partly crystalline feature, consisting of vertical, neutral 
planes, inside which the anions and cations arrange in 
a roughly square lattice, albeit with totally disordered 
tail orientations m- Further decrease of D caused this 
nearly solid IL to squeeze out by successive pairs of ion 
planar layers, positive and negative, thus preserving lo¬ 
cal charge neutrality, reducing the extension of vertical 
planes. 

Starting with large D and an initially fluid IL between 
the plates, the force < F{D) > resisting squeezout, ini¬ 
tially zero, rises as D drops. The IL is partly squeezed 
out, and the part that remains trapped acquires a struc¬ 
ture consisting of Niayer layers (odd for equally charged 
plates, even for opposite charging). The squeezout transi¬ 
tions Niayer Niayer ~ 2 are dear first order transitions 
in this picture, the free enthalpy jumping across barri¬ 
ers between successive H{D) local minima. Decreasing 
D, and thus increasing < F{D) > the lowest enthalpy 
state moves to a lower and lower number of layers. Each 
of these states with progressively lower, integer number 
of IL layers provides a starting state for our successive 
sliding simulations. 

Upon checking, we actually found that configurations 
generated by fast squeezouts often possess internal layers 
that are only partly populated. This kind of ill-prepared 
states can survive as metastable for a long time. When 
faced with these ill-prepared states, initial particle config¬ 
urations were further refined until a sufficiently reliable 
equilibrium particle population was reached. That al¬ 
lowed us to discard as a rule metastable initial configura¬ 
tions, and to describe the sliding behaviour of the stable, 
history independent ones, striving to discard metastable 
states as much as possible within our modest simulation 
durations. 


SIMULATED SLIDING 

As mentioned earlier, there are two types of plate 
charging relevant for sliding friction. The first type, and 
the one generally realized experimentally, is single-plate 
charging relative to an electrode in electrical contact with 
the IL. In this case the aspect that matters is the IL 
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FIG. 2: (Color online) Interplate force F[z) as a function of 
distance, and free enthalpy H(D) calculated by Eq.(l) in the 
absence of sliding and for indicated values of normal load Fn 
and an equal charge on both plates of g = —A^CIcm?. 
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structuring close to the charged plate, whereas the second 
plate, or tip, often insulating, acts as a merely mechan¬ 
ical probe of that structure. We replace this situation 
by simulating two equally charged plates, neutralized by 
adding the right amount of extra ions to the IL. The sec¬ 
ond type, whose experimental feasibility we ignore, but 
which looks even more interesting, is that of opposite 
plate charging. This we simulate without difficulty, re¬ 
stricting the charge magnitude to values small enough to 
correspond to an interplate voltage below 5 V. 

We start off our sliding simulations with equally 
charged plates, endowed with a negative charge density 
q = — Yl\xCjcrv? (about 1/3 of the nominal value quoted 
for mica in experimental conditions) and with the IL con¬ 
fined between them in a variety of Niayer configurations, 
obtained as a result of squeezout down to different D 
values, as detailed in the previous section, further equi¬ 
librated under an applied load Fjq, chosen to be close to 
< F{D) >. Sliding friction depends upon the IL layer 
number Niayer and the load F^. 

Figj^ illustrates two sliding simulations starting from 
Niayer='^ with either Fn = 0.2 nN and Fn = 0.9 nN. 
The simulated system at Fn = 0.2 nN maintains 7 layers 
confined in between the plates, while that at Fn = 0.9 nN 
undergoes a 7 —>■ 5 relayering transition, whose kinetics 
seems favored by sliding. 

As shown in FigHjD, this transition corresponds to the 
sudden drop of the gap D, associated with the expulsion 
of the pair of layers, replacing the Niayer=7 metastable 
state with the stable Niayer=^- This is an example of 
how the lateral driving can facilitate the squeezout of 
layers by mechanically perturbing the IL. 

The close proximity of the 7 —>■ 5 transition can be pre¬ 
dicted by free enthalpy curves, displayed in Fig|^. At 
Fn = 0.9 nN the energy barrier separating the states 
with 7 and 5 layers drops. The role of sliding appears 


FIG. 3: Time evolution of a) the frictional force < Fl > and 

b) interplate gap width D as a function of time for two normal 
loads Fn = 0.2 nN (stable) and Fn = 0.9 nN (metastable). 

c) Free enthalpy curves from squeezout before sliding for in¬ 
dicated values of Fn- Note a stronger local minimum at 
Niayer = 5 for Fn = 0.9 nN, justifying the 7 —>■ 5 transi¬ 
tion under sliding. 

to be kinetic, facilitating the jump to the more stable 5 
layer state. 

Niayer- DEPENDENT FRICTIONAL SHEAR 
STRESS 

The calculated force-distance curves of Ref. [TS] show 
that a single normal load Fn generally gives rise to 
several long lived states with a different number Niayer- 
For example Fn = 0.5 nN can support Niayer = 3, 5, 7 
m, very much as it happens in experiment mini HI]. 
Figure shows the lateral force < Fl > obtained 
as a function of time for different Niayer values, with 
Fn = 0.5 nN and q = —Yl\xClcrv?- Static free enthalpy 
curve in Fig|^ demonstrates that Niayer = 7 is stable 
configuration while 3 and 5 are metastable. A strong 
friction increase is observed with decreasing Niayer, 
in agreement with expectations and with SFA data 
EIIIIS]. The detailed load dependence of the frictional 
shear stress obtained for each fixed number of layers 
is more interesting, and is illustrated in Figj^. Here 
the growth of shear stress with increasing load Fn is 
surprisingly modest. While that might appear at first 
sight in disagreement with experiment, which show a 
nearly linear increase of kinetic friction with load |21) . 
we must recall that the contact area is not a constant 
in SFA, whereas shear stresses require normalizing to 
the area A- Assuming for SFA two perpendicularly 
oriented elastic cylinders in Hertzian contact, the 
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FIG. 4: a) Lateral force < > as a function of time for plate 

charge density of g = —12/cm^ and for different number 
of confined layers Niayer under the same load Fn = 0.5. b) 
Static free enthalpy curve for Fn = 0.5 nN, showing that 
Niayer = 7 is Stable configuration while 3 and 5 are metastable 
ones. 


Numerical simulations Experimental data 



FIG. 5: a) Average shear frictional stress < a >=< FlIA > 
as a function of load Fn for the indicated number of confined 
layers and a charge density of g = —12^C/cm^ on both plates. 
The colored markers refer to the force profiles in panel a) of 
Figi The cartoons suggest how, upon reducing the number 
of layers, the shear plane moves closer to interfaces, where 
particles are strongly bound to the plates, b) Approximate 
behaviour of the shear stress observed with [CdClPyrr] [NTf2] 
confined between two mica sheets in a Surface Force Balance 
experiment. This plot has been obtained by rescaling kinetic 
friction force data from Ref. [21] by a contact area taken as 
A oc Ft^N (ge 0 text). Theoretical and experimental shear 
stresses agree remarkably well, showing a relatively modest 
increase with load. 
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FIG. 6 : Gharge density g = — VljiCIcvn? on both plates, top 
plate driven with velocity V = 2.2m/s. Layer-by-layer den¬ 
sity profiles ( (a), (c) and (e)) and sliding velocity ( (b), (d) 
and (f) ) relative to the bottom plate, for Niayer=7, Niayer=5 
and Niayer=S respectively. Red, blue and green curves repre¬ 
sents cation, anion and tail densities respectively. Gompari- 
son with average expected velocities (dashed lines) shows that 
shear is maximal in the three inner layers (the ’’shear band”), 
moving closer to the plates as Niayer decreases. Error bars in 
panels b, d, f have been reduced ten times for clarity. 


area A oc where Ft = Fadhesion + Fn is the 

effective force between the two sliders. Extracting an 
approximate adhesion force by linearly extrapolating 
the experimental kinetic friction to zero in Fig.l of 
Ref. [^ , we recover the ’’experimentally derived” shear 
stress vs load curves in Fig [^. Remarkably, they now 
show a weaker experimental dependence on load, actu¬ 
ally very close to that predicted by simulations in Fig. [5}i. 

To understand the strong dependence of < Ft > on 
the number of layers illustrated in Fig|^, we now take 
a closer look to the lubricant shear velocity distribution 
during sliding. Fig|^ shows the layer-by-layer resolved 
ion density and velocity profiles along 2 in the liquid for 
7, 5 and 3 layers, Fn = 0.5 nN and V = 2.2 m/s in 
correspondence of green, blue and red symbols in Fig[^ 
respectively. For Niayer = 7 the velocity profile devi¬ 
ates from the linear behavior (red dashed line) indicat¬ 
ing that the layers at the boundaries are bound to the 
plates and tend to move at their same velocity. The plate 
charging promotes the IL in immediate contact to nearly 
solid or at least to a much higher viscosity. Moreover, 
for this relatively high level of charging, the density ra¬ 
tio of first-layer cations over the anode, and of anions 
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over the cathode, is about 1.4, that is only slightly larger 
than one. [14] Both elements support a good grip by 
the plates. The three inner layers constitute by contrast 
the more liquid, effective shear band where the sliding 
is concentrated, and where the coalescence between the 
two solid-like structures protruding from opposite plates 
is easiest to fracture. By reducing the number of layers to 
5 and then to 3, the shear band moves closer and closer 
to the boundaries and even the more solid-like structures 
closer to the plates must break up and move, as shown 
in Fig. ii, f. As a consequence the lateral force < Fl > 
increases. A pictorial representation of this behavior is 
provided by the cartoons of Figj^ Surprising at first 
sight is the seemingly ’’viscous” friction for Niayer = 3, 
with a linear speed profile distribution between the three 
layers, the absence of strong stick-slip as seen in in Figj^ 
and, as we shall see later in Figj^, a linear friction in¬ 
crease with speed. That behavior coexists with a very 
large shear stress, and a very solid structure of this ex¬ 
tremely thin lubricant film. The evolution of layer-by- 
layer relative velocities shown in Fig.6 indicates that the 
yield stresses of the plate/IL interface and of the mid-film 
intra-IL interface, which for Niayer=7,5 were quite differ¬ 
ent, have now become comparable for Niayer=^- As a re¬ 
sult all interfaces slide by about the same amount. The 
meandering wall {x, y) static structure adopted by the 
IL in this pressured-and-charged regime m now inces¬ 
santly ruptures and rearranges during sliding, leading to 
a z-uniform velocity gradient, reminscent of fluid lubrica¬ 
tion, despite its nature very different from that of a fluid. 
The picture changes drastically by reducing the charge 
to <7 = —AfiC/crri^ symmetrically on both plates. The 
mid-film shear band disappears and the slippage occurs 
mostly at the interfaces between the IL and the plates, 
as shown in Figj^ While for q = —12/iC'/cTO^ the top 
plate moved by rupturing the solid-like structures reach¬ 
ing across plates, at this lower charge the film behaves 
effectively as a solid block confined between two slippery 
surfaces. The shearing concentrates mostly at the two 
interfaces, whence the shear stress depends much more 
on the plate-IL interaction. 

The frictional shear stress < a >=< > /A is about 

an order of magnitude smaller than for higher charge and 
displays a much weaker dependence on Niayer ■, as seen 
by comparing Figj^ with Figj^. The rise of average 
frictional force < Fl > for Niayer = 3, displayed in Fig|^ 
(red line) is associated with a stick-slip regime, as shown 
by instantaneous spring forces in Fig[^ corresponding 
to the points indicated in Figj^. This is in contrast to 
the case of high charge, where friction was large but the 
sliding was smooth. 

The essence for this difference of behavior is explained 
by the presence, at low but not at large plate charging, of 
overscreening, a phenomenon well demonstrated in ear¬ 
lier work M- At the small charging of |q| = AyLC/cm^ 
we find that the charge density ratio of first-layer cations 
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FIG. 7: Charge density q = —AyCIcm'^ on both plates, top 
plate driven with velocity V = 2.2m/s. Layer-by-layer den¬ 
sity profiles ( (a), (c), (e) ) and sliding velocity ( (b), (d), 
(f) ) relative to the bottom plate, for Niayer=7, Niayer=5 
and Niayer=S respectively. Red, blue and green curves repre¬ 
sent cation, anion and tail densities respectively. For this low 
charge density the sliding occurs mostly at interface between 
the IL film and the plates. Error bars in panels b, d, f have 
been reduced ten times for clarity. 

over that of the plate anode has in fact risen to about 2.3, 
much larger than one. The overscreened plate is covered 
by a dense counterion blanket, making the plate-IL adhe¬ 
sive grip quite effective. The IL film, on the other hand, 
is solid and even harder to fracture. The plate-IL inter¬ 
face yields under sliding and that causes stick-slip. 


VELOCITY DEPENDENCE OF FRICTION 

After the charge, load, and layer number dependence of 
confined IL friction, we close the symmetric plate charg¬ 
ing part of this study by investigating the velocity depen¬ 
dence. A starting consideration here is that experimental 
sliding velocities in SFA and AFM are very small, typi¬ 
cally reaching 10“®m/s, a value which is at least 6 orders 
of magnitude lower than velocities realistically accessi¬ 
ble in atomic and molecular level simulations. This is a 
standard difficulty, well known and amply discussed [31] ; 
one can in fact learn enough even by simulating at very 
high speeds. Essentially, friction can be either smooth 
and viscous, in which case the shear stress is small and 
proportional to velocity, or it can be stick-slip, in which 
case the shear stress is large and very poorly dependent 
on velocity. In both cases, very crudely speaking, even 
a fast simulation can convey the necessary information 
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FIG. 8: a) Frictional force Fl as a function of time for equal 
and low plate charge density q = —AfiC/crn^, and for different 
number of confined layers Niayer- b) Average shear stress 
a =< Fl > /A as a. function of Fn for indicated number 
of confined layers and A = 177nm^. Colored markers refer 
to the time profiles shown in panel a). The cartoon shows 
that for this low charge density the slippage mostly occurs at 
interfaces between the liquid and plates. 
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velocity independent shear stress; c) quantitative aspects 
are controlled by parameters such as the slider’s masses, 
the pulling spring constant, etc., whose assumedvalues 
are largely arbitrary. 

With these provisions we investigate the dependence 
of IL friction upon driving velocity V, at fixed normal 
load Ftv = 0.5nN. Figj^ and [T0| show the frictional shear 
stress as a function of V for equal-sign charge densities 
of q = —\2^C/cw? and q = —A^C/crn?' respectively and 
for Niayer = 3,5,7. All results for q = —Y2.\xClcm} in 
Figj^show a linear increase of shear stress with velocity 
indicating a viscous-type friction. As suggested earlier, 
this is most likely attributable to overscreening of both 
plates by the extra anions. However < cr > vs F does 
not extrapolate to 0, so that even at F ~ 0 a finite force 
is needed to set the system in motion. This is the force 
needed in order to rupture the vertical IL ’’solid” walls 
that straddle across the two charged plates (Fig[^. 

Also for low charge density q = —d/iC/cm^, the shear 
stress for Niayer = 5,7 grows proportional to velocity but 
here it extrapolates to zero as F —>■ 0 (see FigpTjo). This 
is the hallmark of true viscous friction, and it agrees with 
the top plate moving in a clear regime of smooth sliding. 
The sliding regime drastically changes for Niayer = 3 and 
low velocity F = 1.1m/s, where a clear stick-slip regime 
takes over. The alternation of sticking periods followed 
by rapid and highly dissipative motion of the top plate 
is responsible for the strong increase of average lateral 
force at low velocities, as demonstrated in FigjTO^. By 
increasing F, the stick-slip first turns chaotic and even¬ 
tually disappears at F = 5.5m/s, where smooth sliding 
is recovered. Particle trajectories clearly show this tran¬ 
sition. At low F the IL structure retains its worm-like 
meandering planes or chains, which we described in ref. 
m, and which were recently observed in experiments, 
m and behaves essentially as a solid confined between 
two slippery planes. 

The transition from stick-slip to smooth sliding is 
shown in Fig |10^ , reporting two force profiles for 
Niayer = 3 Corresponding to driving velocities indicated 
by arrows in Fig|10|3. 


FIG. 9: a) Frictional force Fl as a function of time for 
Niayer = 5, q = —12)10/071^ and velocities indicated by ar¬ 
rows in panel b. b) Average shear stress < a >=< Fl > /A 
as a function of driving velocity V for different number of con¬ 
fined layers Niayer = 3, 5, 7. The value of normal load here is 
Fjv = 0.5nN. 

(the viscous friction coefficient in the first, the full shear 
stress in the second). Of course it will always be neces¬ 
sary to bear in mind that a) increasing velocity may turn 
stick-slip into smooth sliding; and b) microscopic scale 
stick-slip does not necessarily show up as macroscopic or 
mesoscopic oscillations, and is generally detectable as a 


FRICTION WITH OPPOSITE PLATE 
CHARGING 

In this final part we study how friction force is affected 
by opposite charging of the two confining plates, such 
as one would obtain by applying a DC voltage between 
them. Probably not that easy to realize experimentally, 
this configuration is nonetheless quite interesting. Upon 
balancing the opposite charge on the two plates, the equi¬ 
librium layering structure now forms an even number of 
layers, a rearrangement which preserves local charge neu¬ 
trality. In order to explore expeditely the dependence 
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FIG. 10: a) Lateral force Fl as a function of time for 3 layers, 
q = —AfiCIcm^ and velocities indicated by arrows in panel 
b. b) Average shear frictional stress < a >=< Fl > IA as a 
function of driving velocity V for different number of confined 
layers Niayer- The value of normal load here is Fn = 0.5nN. 

of friction upon the magnitude of the charge imbalance, 
we conduct a time-dependent simulation where we pe¬ 
riodically, but very slowly and as adiabatically as pos¬ 
sible, modulate in time the on-plate charge in the form 
g = ±Qsin(27r</r), where the two opposite signs refer 
to the two opposite plates. Adiabaticity is pursued by 
choosing a period r large enough compared to the char¬ 
acteristic time scale of all relevant ionic rearrangements. 

The magnitude of charging used in this work, Q < 
l2^C/cm? is dictated, as discussed previously [T5] . by 
the requirement that, assuming an average dielectric con¬ 
stant of 2, the plate-plate voltage should not exceed 5V, 
taken as a reasonable upper limit in an experimental sit¬ 
uation. In this respect therefore, our study does not over¬ 
lap with previous work, where in a different model |13j 
the behavior at Q > 16^C/cw? was investigated. 

Assuming a constant load Fpf = 1 nN, a charge oscilla¬ 
tion period t = 12 ns, a magnitude Q = jcrr? and 
a sliding velocity V = 22m,js we obtain the spontaneous 
evolution of the system shown in Figj^ What happens 
is, in a nutshell, that the interplate distance D sponta¬ 
neously changes, and so does the sliding friction Fl as 
the charge changes. The simulation begins with a charge 
density of |g(t = 0)| =0 on plates and the IL arranges in 
one single neutral layer indicated as In in Fig |ll[ which, 
as shown in the static squeezout enthalpy of Fig{^ is the 
most stable configuration under the fixed applied loading 
force Fn = 1 nN and low charge. As the plate charge 
increases, the liquid in the lateral reservoirs is sucked in 
the gap (despite a constant load), reaching a four layer 
solid state. Actually, we expected it to go back to six 
layers, which is the equilibrium free enthalpy minimum, 
as shown in Fig|12[ It probably would if we could simu- 



time(ns) 


FIG. 11: a) Frictional sliding evolution while slowly changing 
opposite charges on the two sliding plates. Parameters: Initial 
Niayer = 1, Fn = 1 uN, charges q = ±Q sin(27rt/T) with 
Q = lOnCIcnF and r = 12 ns, sliding velocity V = 22mls. a) 
Spontaneous evolution of plate-plate distance D. b) Number 
of particles Npart in the interplate gap. c) Sliding force Fl 
as a function of time, d) Time evolution of top and bottom 
plate charge. Gharging causes the plate separation D to open 
up at charge maxima (IL suck-in) and to collapse at charge 
zeroes (IL squeezout), all under constant load. The inset in 
panel a) shows a side view of the confined liquid structure for 
4 layers {q = Q) and for the single neutral planar layer In 
{q = 0). Friction is large when |g| and voltage are maxima, 
and charging solidifies the IL inside the gap D. Gonversely, 
when q is near zero, the IL melts, the plates collapse together 
squeezing out all fluid except for a neutral monolayer where 
both ions pack in a square lattice, resulting in a very low 
friction. 



FIG. 12: (Golor online) Enthalpy curves for indicated values 
of charge density on plates. Numbers indicated in figure de¬ 
note metastable states at different number of layers. Peak 
heights decrease with charge and arrows indicate charge in¬ 
duced transitions between the states In and 4. 
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late for much longer times, but it does not on our limited 
time scale, due to tha strongly reduced mobility of ions 
in the near-solid film. This strong and rigid four layers 
structure softens as the plate charge decreases, eventu¬ 
ally turning more fluid and squeezing out reaching again 
the single neutral layer arrangement In. These charging- 
induced relayering and squeezout are accompanied by 
drastic change in frictional behaviour of the liquid film. 
The large shear stress at high plate charge (high voltage, 
near 5 eV) is replaced by a nearly vanishing shear stress 
at zero charge. Here therefore, contrary to the equal- 
charge case of the previous sections, the thinnest film 
yields the lowest friction, rather than the opposite. This 
is because, independently of sliding, the plate charging 
induces structural transformations of freezing/melting of 
the the IL film|15|. The frozen IL supports the applied 
load for large charge but the liquid IL does not for zero 
charge. The IL film structure at g ~ 0 shrinks to a sin¬ 
gle, planar, square and neutral IL monolayer. Here the 
neutral tails stick out of the plane, effectively lubricat¬ 
ing the sliding. A qualitatively similar mechanism where 
alkyl tails stick out of a well-packed anion-cation plane 
has been proposed in a recent experimental study [S] to 
explain the potential-dependent frictional behaviour of 
an AFM tip in presence of an IL. 

We did change and explore a variety of parameters, in¬ 
cluding charging magnitude, applied load, and sliding 
speed in order to test the main aspects of these results. 
The processes of melting and freezing, and of expulsion 
and sucking always remained very asymmetrical between 
the fast squeezout, associated with sudden melting of the 
IL film upon decreasing jgl towards zero, and the more 
gradual sucking-in and freezing, attained upon increas¬ 
ing l^l from zero in FigpT^. The electrically driven fast 
emptying but slow refilling of the interplate gap could 
be used, together with some spatial asymmetry to build 
an effective effective electro-pumping of ions in and out 
of the gap. The frictional shear stress shows a very 
large charge dependence, as a direct consequence of the 
changes of structure and of gap width D{q) caused by 
charging. Although noisy, the large shear stress at charge 
maxima drops by an order of magnitude or more follow¬ 
ing squeezout at zero charge. This strong effect appears 
promising in view of experimental verifications, and of a 
possible use for the control of friction. 

Lastly, we briefly dwell on questions connected with 
work, heat and temperature, which come about with re¬ 
spect both to charging-induced melting and freezing, and 
to frictional work and heating. In the absence of sliding, 
electrically induced freezing and melting would be asso¬ 
ciated with the absorption or release of latent heat. That 
would be large if the solid was ordered, e.g., crystalline 
or nearly crystalline, as is the case for example in the 
NaCl model of the IL m- The more the solid resem¬ 
bles a glass, however, the smaller the latent heat will be, 
owing to structural near-identity between glass and liq- 



time(ns) 

FIG. 13: (Color online) Comparison between the charge in¬ 
duced squeezout and relayering processes in the presence and 
absence of sliding. Spontaneous evolution of plate-plate dis¬ 
tance D ((a) and (c)) and temperature ((b) and (d)) in com¬ 
parison with time evolution of top and bottom plate charge 
(e). In the absence of driving the temperature remains stable 
at the set value of T = 225K, while in the presence of driving 
it temporarily deviates from this value. 

uid. That is the case of our TM model, where indeed the 
latent heat of melting was very small in the simulated 
caloric curves m- Without sliding, the presence of a 
latent heat is undetectable in electrically induced freez¬ 
ing and melting cycles shown in FigfT^,d. Of course, 
these data were obtained with a thermostat, canceling 
any large temperature fluctuations: despite the themo- 
stat, dips and peaks would be well visible in presence of 
a robust, crystal-like latent heat. Frictional heating is 
instead quite visible in the sliding results under cycled 
charging of FigfT^,b. Despite the thermostat, the tem¬ 
perature is seen to oscillate, downwards when the solid 
confined IL lubricant goes liquid at g = 0, up when it 
goes back to solid at |q| = maximum and where there¬ 
fore the frictional Joule heating is much stronger, down 
again upon melting where the friction drops once again. 
Because of the presence of the thermostat we did not at¬ 
tempt a quantitative connection between the frictional 
work and the residual temperature oscillation, but the 
connection between the two is quite clear. 


CONCLUSIONS 

The electrical charging of plates in nearly atomic con¬ 
tact influences an IL when trapped in a nanoscale sized 
gap between them. As a result, the frictional shear stress 
realized upon mutual sliding of the plates is modified by 
the charging. The present simulation study conducted 
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within the simple ’’tailed model” previously developed 
for the IL explores some of the rich variety of phenom¬ 
ena that can take place in this context. The IL becomes 
structured into alternating charge layers, whose number 
is odd when the plates are equally charged, even when 
they are oppositely charged. The trapped IL film may 
develop a solid-like rigidity, with different and opposite 
frictional consequences depending on charging level and 
other parameters. When the plate grips the solid-like film 
and the shearing takes place in its middle the friction is 
large; but when the grip is less effective the shear con¬ 
centrates at the plate-film interface and friction drops. 
Alternatively, at zero charging and under the same con¬ 
ditions the ionic liquid can melt, thus lubricating and 
dramatically reducing the sliding friction. 

Even though the parameter-dependent nature of our 
results does not make a comparison with any specific 
experimental system particularly compelling, we believe 
that several elements uncovered either connect with ob¬ 
servations already available, or else suggest pursuing 
newer ones. The modest dependence of friction upon 
load - unusual in the context of classical friction - is 
among the former. The possibility that electrical charg¬ 
ing, in whichever manners that could be realized, could 
influence the tribology of confined IL by causing its elec¬ 
trically driven freezing and melting is among the latter, 
and deserves to be tested, supporting the possibility of 
an electric control of IL lubricated friction. 
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